home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
CD ROM Paradise Collection 4
/
CD ROM Paradise Collection 4 1995 Nov.iso
/
program
/
nrpas13.zip
/
CHSTWO.DEM
< prev
next >
Wrap
Text File
|
1991-04-29
|
1KB
|
54 lines
PROGRAM d13r9(input,output);
(* driver for routine CHSTWO *)
CONST
nbins=10;
npts=2000;
TYPE
barray = ARRAY [1..nbins] OF real;
VAR
glinext,glinextp : integer;
glma : ARRAY [1..55] OF real;
chsq,df,prob,x : real;
i,ibin,idum,j : integer;
bins1,bins2 : barray;
(*$I MODFILE.PAS *)
(*$I GAMMLN.PAS *)
(*$I GSER.PAS *)
(*$I GCF.PAS *)
(*$I GAMMQ.PAS *)
(*$I RAN3.PAS *)
(*$I EXPDEV.PAS *)
(*$I CHSTWO.PAS *)
BEGIN
idum := -18;
FOR j := 1 to nbins DO BEGIN
bins1[j] := 0.0;
bins2[j] := 0.0
END;
FOR i := 1 to npts DO BEGIN
x := expdev(idum);
ibin := trunc(x*nbins/3.0)+1;
IF (ibin <= nbins) THEN bins1[ibin] := bins1[ibin]+1.0;
x := expdev(idum);
ibin := trunc(x*nbins/3.0)+1;
IF (ibin <= nbins) THEN bins2[ibin] := bins2[ibin]+1.0
END;
chstwo(bins1,bins2,nbins,-1,df,chsq,prob);
writeln;
writeln('dataset 1':15,'dataset 2':15);
FOR i := 1 to nbins DO BEGIN
writeln(bins1[i]:13:2,bins2[i]:15:2)
END;
writeln;
writeln('chi-squared:':18,chsq:12:4);
writeln('probability:':18,prob:12:4)
END.